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We extend the numerical renormalization-group method to treat Bose-Fermi Kondo models 
(BFKMs) describing a local moment coupled both to a conduction band and to a dissipative bosonic 
bath representing, e.g., lattice or spin collective excitations of the environment. We apply the method 
to the Ising-symmetry BFKM with a structureless band and a bath spectral function rj(uj) oc uj s . The 
method is valid for all bath exponents s and all temperatures T. For < s < 1, the range of interest 
in the context of heavy-fermion quantum criticality, an interacting critical point, characterized by hy- 
perscaling of exponents and w/T-scaling, describes a continuous quantum phase transition between 
Kondo-screened and localized phases. For Ohmic dissipation s = 1, where the model is relevant to 
certain dissipative mesoscopic qubit devices, the transition is found to be Kosterlitz-Thouless-like. 
In both regimes the impurity spectral function for the corresponding Anderson model shows clearly 
the collapse of the Kondo resonance at the transition. Connection is made to other recent results 
for the BFKM and the spin-boson model. 
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I. INTRODUCTION 

Quantum impurity models that incorporate dissipative 
effects of the environmenfeiiS are highly topical. Applica- 
tions of such models, which represent the environment 
by one or more bosonic baths, are numerous and in- 
clude noisy quantum dotSfiii decoherence of qubits in 
quantum computation^S^^ and charge transfer in bio- 
logical donor-acceptor systemsii2iii Furthermore, within 
dynamical mean-field theory (DMFT) and its exten- 
sions, certain lattice models, of interest, e.g., in rela- 
tion to heavy-fermion quantum criticality or the inter- 
play of electrons and phonons in strongly correlated elec- 
tron systems, also reduce to impurity problems of this 
nature^ 2 *^ In this paper we provide a comprehensive 
account of the Ising-symmetry Bose-Fermi Kondo model 
(described further below) obtained via an extension of 
Wilson's numerical renormalization-group method^ 

The simplest example of a quantum impurity coupling 
to a bosonic bath is the celebrated spin-boson model 
(SBM), describing the dynamics of a two-state system 
coupled to a bath of harmonic oscillators. For a review 
see Ref. QJ. The essential physics is the competition be- 
tween the amplitude for tunelling between the two spin 
states S z — ±| (leading to a "delocalized phase") and 
the effects of the environment embodied in the coupling, 
via S z , to the bosonic bath, which tends to "localize" 
the system in one or other of the spin states. Not sur- 
prisingly, as the generic model for quantum dissipation, 
the SBM itself finds many applications in physics and 
chemistry^ 

More general than the SBM are so-called "Bose-Fermi" 
quantum impurity problems involving both fermionic and 
bosonic baths. The Bose-Fermi Kondo model (BFKM) 
was originally discussed in the context of an extended 
dynamical mean-field theory (EDMFT) treatment of a 
two-band extended Hubbard model 15 and describes a lo- 



cal moment coupled not only to a fermionic conduction 
electron band — the "regular" Kondo model (for a re- 
view, see Ref. 0) — but also to a dissipative bosonic 
bath specified by a spectral function B{oj) that is typi- 
cally modeled as B(uj) cx uo s at low frequencies. 

Of particular interest is the BFKM with a sub-Ohmic 
bath s < 1, which has been studied within a perturba- 
tive (in e = 1 — s) rcnormalization group approach^iiiS 
The coupling of the impurity to the bosonic bath com- 
petes with the Kondo effect and leads, for < s < 1, 
to an unstable interacting critical point that describes a 
continuous quantum phase transition (QPT) between a 
Kondo-screened phase, in which the impurity moment is 
quenched by conduction electrons ((S z ) — 0), and a local- 
ized phase where (S z ) ^ 0. Thus, the BFKM provides an 
example of an impurity QPT, encapsulating the compe- 
tition between local and spatially extended physics. For 
a recent review of impurity QPTs, see Ref. Il9i 

One possible direct realization of the Bose-Fermi 
Kondo model (BFKM) is a noisy metallic quantum dot: 
a large micron-scale dot coupled to a reservoir lead and 
subject to Ohmic dissipation (s = 1) arising via fluctua- 
tions in the gate voltage. Here, the two impurity "spin" 
states represent the two possible charge states of the dot 
close to a degeneracy point. Quantum charge fluctua- 
tions on the dot can be supressed by voltage fluctuations, 
restoring perfect Coulomb-blockade staircase behavior. 3 
The model has also been invoked to study, within a dy- 
namical large- N treatment, the quantum-critical Kondo 
effect in a single-electron transistor coupled to both the 
conduction electrons and spin waves of ferromagnetic 
leads^ 

However, to date the BFKM has received most at- 
tention in the context of heavy-fermion quantum crit- 
icality, arising in an EDMFT treatment of the Kondo 
lattice modeli 21 ' 22 i 23 ' 24 The Kondo lattice model is be- 
lieved to capture the essential physics underlying heavy- 
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fermion QPTs: the dynamical competition between 
Kondo screening of /-shell moments on each lattice site 
by conduction electrons, which leads to a paramagnetic 
heavy-electron metallic phase; and antiferromagnetic or- 
dering of those moments due to the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction. In EDMFT, which 
allows for the nonlocal quantum fluctuations arising from 
RKKY interactions between local moments, all correla- 
tion functions of the Kondo lattice problem can be cal- 
culated from a self-consistent Bose-Fermi Kondo model. 
Here, the bosonic bath embodies the effects, at a given 
site, of the fluctuating magnetic field generated (via the 
RKKY interaction) by local moments at other lattice 
sites. 

The above is of considerable current interest in connec- 
tion with so-called local quantum criticality. In this sce- 
nario of a heavy-fermion QPT, destruction of the Kondo 
effect accompanies divergence of the spatial correlation 
length at the antiferromagnetic ordering transition. Such 
a local quantum-critical point may arise within EDMFT 
- as a critical impurity solution satisfying the self- 
consistency — when magnetic spin fluctuations are two- 
dimensionaliSiiSS The salient non-Fermi liquid features 
of several materials, including CeCus.gAuo.i, which re- 
sist the standard spin-density wave description^ ap- 
pear to be consistent with a critical point of this nature. 
However, whether or not a local quantum-critical point 
(QCP) arises self-consistently under a nonperturbative 
treatment that considers both paramagnetic and antifer- 
romagnetic phases of the model, is an important question 
and one that has not been satisfactorily answered by re- 
cent quantum Monte Carlo studiesi^SLSS These studies 
reach opposing conclusions concerning the nature of the 
T = transition, in part due to the inherent limita- 
tions of accessing the lowest temperature scales within 
quantum Monte Carlo. Similar EDMFT mappings to a 
BFKM have been applied to the disordered Kondo lattice 
model to study interplay of Kondo screening and RKKY 
interactions in electronic Griffiths phases, 2 - and to a t- 
J model of a highly incoherent metal close to a Mott 
transition^ of possible relevance to the pseudgap phase 
of the cuprates. 

Bose-Fermi impurity problems may also arise through 
a conventional DMFT treatment of lattice models that 
include both electron correlation and electron-phonon 
interactions. The Hubbard-Holstein model is the sim- 
plest such model2ii22iS and has been discussed in con- 
nection with cuprate high-T c superconductors, heavy- 
fermion systems, manganites and fulleride superconduc- 
tivity. For a recent discussion of Bose-Fermi impurity 
problems in the context of DMFT, see Ref. H^. 

The preceding remarks serve to motivate the study of 
Bose-Fermi impurity problems as a route to the descrip- 
tion of a variety of physically relevant models of strongly- 
correlated systems. To what extent the impurity models 
(or impurity models plus DMFT self-consistency for lat- 
tice models) capture the relevant physical behavior can 
only be assessed if a reliable impurity solver is employed. 



Such a technique should be capable of handling static 
and dynamical properties for all interaction strengths, 
bosonic bath spectra and temperatures. This is a de- 
manding requirement, even for purely fermionic impu- 
rity models, and many of the standard techniques de- 
veloped for the latter are here simply inapplicable. It 
is natural to turn to Wilson's nonperturbative numerical 
renormalization-group (NRG) method, 14 which has so far 
provided a controlled numerical treatment of the thermo- 
dynamics, dynamical response functions, and transport 
properties of a range of fermionic impurity models and 
of lattice models within DMFT. Recent impurity-model 
applications include studies of quantum criticality in 
the pseudogap Kondo^i and Anderson^ models, uncon- 
ventional physics in double quantum dots,^^^^ and 
non-Fermi liquid behavior in a frustrated three-impurity 
Kondo model4£ Applications to lattice systems include 
studies of the Mott-Hubbard metal-insulator transition 
in the Hubbard model^i and the paramagnetic phase of 
the periodic Anderson model relevant to heavy-fermion 
systems^ For a review, see, e.g., Ref. l43l The scope of 
potential applications of the NRG method has recently 
been widened following a successful treatment of the spin- 
boson model via a pure bosonic NRG^ 4 *^ The latter has 
provided a good account of the critical properties of the 
SBM for both Ohmic and sub-Ohmic bath spectra. 

This paper reports a direct study of the Bose-Fermi 
Kondo model using an NRG method extended to han- 
dle simultaneously both fermionic and bosonic degrees of 
freedom. We demonstrate that our approach provides 
a good account of the critical properties of the Ising- 
symmetry BFKM for both a sub-Ohmic bosonic bath 
with s < 1, as reported briefly in Ref. E^. and for the 
transition in the Ohmic case s = 1. We provide an impor- 
tant confirmation that the Ising-BFKM and SBM belong 
to the same universality class, as implied by bosonization. 
The NRG treatment described is reliable for all bosonic 
bath exponents and temperatures. Importantly, it opens 
the way for a reliable numerical solution of the Kondo 
lattice model by imposing EDMFT self-consistency on 
the BFKM solution, as we have discussed elsewhere^ 

The paper is organized as follows. After the essential 
background to the BFKM is given in Sec.[Hj an outline of 
the NRG treatment of the model — which involves loga- 
rithmic discretization of the energy axis, mapping of the 
Hamiltonian to a chain form, and the iterative numeri- 
cal solution of the problem at a sequence of decreasing 
energy scales — is given in Sec. 11111 

Results arising from the approach are presented in two 
parts, beginning with the case of sub-Ohmic dissipation 
in Sec. IIVI The flow of the NRG transformation and its 
fixed points is discussed in Sec. IIV Al An interacting 
critical fixed point separates a Kondo-screened phase, 
for weak coupling to the bosonic bath, from a strong- 
coupling localized phase with no Kondo resonance. The 
form of the phase boundaries, in particular the relation- 
ship of the critical coupling to the bosonic bath expo- 
nent and the Kondo scale of the pure-fermionic problem, 
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is discussed in Sec. llVBl In Sec. WU for < s < 1, 
we determine the correlation-length exponent that gov- 
erns the vanishing, at the quantum-critical point, of a 
low-energy crossover scale. The response to a local mag- 
netic field applied at the impurity site provides a use- 
ful probe of impurity quantum phase transitions; in Sec. 
IIVDI we discuss the static and dynamical local suscep- 
tibility. Distinct behavior obtains in either stable phase 
and at the QCP. Critical behavior is manifest as power 
laws (with anomalous exponents) in the temperature and 
frequency dependences, which reflect the fact that the 
critical physics is scale-free. We identify an order param- 
eter for the problem, which, for < s < 1, vanishes con- 
tinuously as the transition is approached. Various critical 
exponents are extracted and are shown to obey hyperscal- 
ing relations. The latter, together with w/T-scaling at 
the critical point, are consistent with an interacting fixed 
point. In Sec. HVEl we consider the single-particle spec- 
tral function A(uj), in which the collapse of the Kondo 
resonance at the QCP is directly manifest. The critical 
properties and scaling behavior are discussed in detail. 

The nature of the transition changes qualitatively at 
s — 1, as reflected in the divergence of the correlation 
length exponent as s — > 1~. We consider the Ohmic case 
(,s = I) in Sec. demonstrating that the transition is 
Kosterlitz-Thouless-like, with characteristic vanishing of 
the renormalized Kondo scale and a jump in the order 
parameter at the critical coupling. The spectral function 
and its scaling properties are discussed and compared to 
their sub-Ohmic counterparts. We also briefly examine 
our results in the context of the noisy quantum boxi 3 i 4 i 5 

The paper concludes with a brief summary and outlook 
for future work. 



II. BACKGROUND 

In the remainder of this paper we focus on the Bose- 
Fermi Kondo model with Ising anisotropy in the bosonic 
couplings. This is the case most relevant to the EDMFT 
treatment of CeCug-zAua; and to the noisy quantum dot 
system. We treat the model by extending the NRG 
method to encompass Bose-Fermi impurity problems. 
Though described in Sec. 1 1 1 1 1 for the Ising-BFKM, the 
approach is straightforwardly generalized to the xy and 
isotropic problems. Results for the latter cases will be 
discussed elsewhere. In the following subsection we de- 
fine the model to be studied and identify key quantities, 
before going on to discuss the connection to the spin- 
boson model in Sec. Ill Bl We set h — gps = &b = 1. 



A. Ising-anisotropic Bose-Fermi Kondo model 

The Hamiltonian for the Ising-symmetry Bose-Fermi 
Kondo model is given by 

H = H F + H B +H intl (2.1) 



where 

Hf =J2 €kC la c ^ (2-2) 

k.a 

and 

H B = J2^^i (2-3) 

describe, respectively, the conduction band and a dissi- 
pative bosonic bath. Both are taken to be noninteracting 
and characterized by the densities of states 

p(e) = £<y(e-e fc ) (2.4) 
fc 

and 

*7(w)=5^*(w-w,). (2.5) 

The interaction of the impurity with the baths is con- 
tained in 

H m t = ^JoS- ^2 ej^owcfev +S z ^2g q (<f> q + (j)L q ). 

k,k' .a. a' q 

(2.6) 

The first term describes the interaction of a spin-i local 
moment with the on-site spin of the fermionic band via 
the usual Kondo coupling Jo; the second term describes 
its coupling, of magnitude g q , to the bosonic bath with 
oscillators {4> q }- 

We assume a featureless metallic host conduction band 
of width 2D described by a density of states 

p(e) = p Q e(D - \e\) (2.7) 

with 0(x) the Heaviside step function. The bosonic bath 
is entirely determined by the spectral function 

B(u;)=7rX>g%-^), (2.8) 
q 

which we assume vanishes outside the range < u> < luq. 
Note that the bath density of states rj(cu) and coupling g q 
do not require separate specification. Writing x = e/D 
and y = u>/u>q, the Hamiltonian is conveniently written 
in the following one-dimensional form: 

H F + H B = D dx x c\ a c xcr +ujq dy y 4> y 4> y , 
J-i Jo 

(2.9a) 



Hmt = PppJoS-^ /dgOVtr'/otr' + UJQ^J S z (b + b f Q ) , 
era' 

(2.9b) 

extracting a dimensionless "dissipation strength" Bq 
from B{uj). The dimensionless operators 

1 f +1 

ha = J dx C xa (2-10) 
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and 



6o = 6- 1 / dy W{y)<j> y , 
Jo 



the latter denned by a weighting function 



W(y) 



and a normalization factor 



< B{lj oV ) 

BqLUq 



B 1 



dy W 2 (y), 



(2.11) 



(2.12) 



(2.13) 



respectively annihilate the unique combination of 
fermions and bosons that couples to the impurity. The 
operators c xa and <f> y obey the commutation relations 
{4<7> c *v} = S(x - x')5 aa , and [0t, <j) y ,] = 5{y - y') re- 
spectively, while {fl a , foa'} = <W and [b Q , b ] = 1. 

Of principal interest is the low-temperature physics of 
the model, determined by the form of the bath spectrum 
at low energies. We take B(u>) to be given by 



for < uj < ujq, 
otherwise, 



(2.14) 



with s > — 1 necessarily. The bulk of the paper concerns 
the sub-Ohmic regime s < 1; in particular we focus on 
< s < 1 where the key physics is governed by an inter- 
acting quantum-critical point. However, we also study 
the Ohmic case (s = 1) in Sec.lVl 

To connect with previous work on the BFKM, we will 
assume that g q is g-independent (g q = go) and that the 
bath spectral function is of form 



(Kq/tt) luq s lo s for < u < lu 







otherwise, 



(2.15) 



and hence, from Eqs. Ij2.8|) and Ij2.5jl . B = (K g ) 2 



B. Connection to the spin-boson model 

The spin-boson model is the simplest model of quan- 
tum dissipation, t describing a two-level system coupled 
to a bosonic bath: 



H SB =H B - AS x + eS z + S, 



Em 



+ 0LJ. (2.16) 



Here, A is the amplitude for tunneling between the 
"spin" states S z = ±=, and e represents an additional 
bias. The coupling to the environment is characterized 
by the spectral function J{uj) = Ag<5(cj — uj q ). It 

is well knowri that there exists a correspondence be- 
tween the Ohmic spin-boson model and the anisotropic 
Kondo model (AKM), described by a Hamiltonian H = 



Hp + H int with Hp given by Eqs. J^D, l|2~i)l and (|2~7|l . 

and 



H\ r 



J± 



fc.fc' 



f E( c lT c fe'T - ^i)S z + hS z . (2.17) 



fc.fc' 



Specifically, following bosonization of the AKM (see, e.g., 
Ref-ET 



J{w) 



2Ttauj for co < lo c 



(2.18) 
(2.19) 
(2.20) 
= [1- 



where the dimensionless dissipation strength a 
(2/7r)tan" 1 (7rp J z /4)] 2 and lu c = 2D. 

The Ohmic case of the SBM is well understood. The 
essential physics is a Kosterlitz-Thouless quantum phase 
transition between a delocalized phase for a < a c — in 
which the effect of dissipation is to renormalize the tunel- 
ing amplitude — and a localized phase for a > a c , where 
tunneling between the spin states is absent. For A — > 0, 
a c = 1. Thus, the Kondo scale Tr- of the AKM can 
be interpreted as a renormalization of the bare tunnel- 
ing amplitude poJ± due to frictional effects of the elec- 
tronic environment (particle-hole excitations) controlled 
by poJz- The antiferromagnetic phase of the AKM corre- 
sponds to the delocalized phase of the SBM, while the fer- 
romagnetic sector corresponds to the localized phase. For 
the regular isotropic Kondo model, J z = J± = Jo, the 
corresponding SBM dissipation strength takes the value 
a — 1~ for po J <C 1. 

The SBM with sub-Ohmic dissipation has only recently 
been studied in detail, using perturbative* and numer- 
ical renormalization group^ methods. For bath expo- 
nents < s < 1 these studies establish the existence of 
an interacting quantum-critical point describing a con- 
tinuous QPT between delocalized and localized phases, 
with a c oc A 1 ^ for A <C lu c . 

Recent studies of the Ising-symmetry Bose-Fermi 
Kondo model with Ohmic dissipation^ have exploited 
the connection between the AKM and Ohmic SBM. For 
example, the essential results of Ref. ^5| have been dis- 
cussed in the context of the Ohmic BFKM, 4 of relevance 
to certain dissipative mesoscopic qubit devices. More 
generally, the Ising-symmetry BFKM can be mapped to 
a spin-boson model whose spectral function J(uS) com- 
prises an Ohmic contribution from the bosonized con- 
duction band and a contribution B(u>) oc lj s from the 
BFKM bath. The resulting J(uo) is dominated by B{ui) 
at low frequencies and thus it is to be expected that the 
Ising symmetry BFKM and SBM have identical univer- 
sal critical properties for s < 1. However, the presence 
of the Ohmic contribution acts to renormalize the bare 
tunneling amplitude entering the model, such that the 
sub-Ohmic BFKM corresponds to a SBM, Eq. 
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with J(oj) oc uj s for uo <C lu c and the bare Kondo scale 
Tk (not J) the appropriate tunneling amplitude. This is 
borne out by the numerical results of Sec. IIVI 



III. NRG TREATMENT 



In Wilson's original NRG treatment of the Kondo 
model 14 - [as described by H = Hp + H int , Eqs. J23, 
(|2.6|) with g q = 0, and 1|2.7[) ] the conduction band con- 
tinuum of energies e is sampled at the discrete set of 
energies e = ±Z?A~ m for m = 0,1,2,..., where A > 1 
parametrizes the discretization. By a mapping of Hp 
to a tight-binding Hamiltonian using the Lanczos pro- 
cedure, the Kondo Hamiltonian is cast in the form of a 
semi-infinite fermionic chain with the impurity coupled 
only to the end site n = 0. Owing to the discretization, 
the tight-binding coefficients decay with site number as 
DA~ n for large n, which feature allows for an iterative 
solution to the problem by numerically diagonalizing the 
Hamiltonian for progressively longer finite chains. The 
procedure amounts to including excitations on an expo- 
nentially smaller energy scale at each step, with the full 
Hamiltonian being recovered as the limit 



H= lim ah~ N/2 DH 



N- 



N- 



(3.1) 



Here, Hm describes an (N + l)-site (fermionic) chain, 
rescaled so that the lowest energy is of order unity and a 
meaningful comparison between the {Hn} can be made. 
The constant a is conventional and given by a = ^ Ai (1+ 
A" 1 ). 

The renormalization-group transformation R, defined 
by Hjy = RHn-i, eventually reaches a scale-invariant 
fixed point of R that determines the low-temperature 
properties of the system. At a fixed point, the low- lying 
eigenstates and matrix elements of any observable oper- 
ator are identical for iterations N and N — 2. That the 
fixed points are of the transformation R 2 and not R re- 
flects the fundamental inequivalence of fermionic chains 
containing odd and even numbers of sites, as discussed, 
e.g., in Ref.l49l 

In the following we outline the transformation of the 
BFKM, Eqs. and (2.9), to the form required for 

an NRG treatment. Including a coupling between the 
impurity and a bosonic bath introduces two additional 
considerations. First, and most obviously, the bosonic 
Hilbert space must be truncated from the outset, intro- 
ducing a bosonic truncation parameter. The effects of 
this truncation on the critical properties of interest must 
be systematically studied. Second, the iterative proce- 
dure employed must preserve the spirit of the NRG ap- 
proach, taking into account all excitations of the same 
energy scale at the same iteration. We turn to a discus- 
sion of these issues in the following section. 



A. Discretization and mapping 

The bosonic bath [conduction band] continuum of en- 
ergies < w < ujq [— D < e < D] is replaced by a discrete 
set: woA~ m [±DA~ m ] for integer m > 0. The main 
approximation of the NRG is to select just one bosonic 
[fermionic] state from each energy interval. We retain 
the single state that couples directly to the impurity such 
that bo in Eq. (|2.11|) is given exactly by 



6o = B" 1 ^B, 



(3.2) 



with 



and 



K = S" 1 / dy W(y) 4> y (3.3) 

j A _ (m+1) 



Bl = / dy W 2 (y). (3.4) 



For the fermions, 



V <m 



(3.5) 



The operators c^ na and c< CT arc defined for the mth pos- 
itive and negative energy intervals respectively and are 
given by 



dx 



A -Cm+1) 



dx C X a- 



(3.6) 
(3.7) 



The terms Hb and Hp are then approximately written 



as 



k>m 4 > ln4 > r. 



(3.8) 



Hp ~ D (e^ C^tr c ma + e m C ma C ma)i (3-9) 



where 



w m — B m 



A-(m+l) 



dy W(y)y, (3.10) 



= ±i (A- m + A~( m+1 )) , (3.11) 



and the approximation is seen to amount to selecting 
the state of weighted mean energy from each interval. 
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fermions 



fo h : h h U 




bosons 



H F = [tnfLfna + T n (f^ a f n -l <<7 + fl_ 1>(7 fna) ■ 

(3.15) 

The tight-binding coefficients e n , r n , e„ and t n (n — 
0,1,2,...) encode all information about the conduction 
band and bosonic bath and are determined numerically, 
for given B{ui) and piui), via Lanczos recursion relations. 
For example, e ra and t n follow from numerically iterating 



FIG. 1: Schematic of the chain form of the BFKM Hamilto- 
nian. Fermionic sites of the tight-binding chain are shown 
as small circles, while bosonic sites are shown as squares. 
Only the end site of each chain couples directly to the im- 
purity spin. The coupling between the fermionic sites decays 
as T n ~ A -n//2 for large n, while the bosonic couplings fall off 
as t n ~ A~™. The dashed boxes enclose, from innermost to 
outermost, the sites included in iterations 0, 1, and 2, respec- 
tively. 



The infinite number of states neglected in each inter- 
val are orthogonal to the state retained and couple only 
indirectly to the impurity, which coupling vanishes as 
A — > 1. As is standard in the NRG approach, we ne- 
glect this coupling for A > 1 and drop the irrelevant 
decoupled constant term that thereby results. The dis- 
cretization error incurred vanishes as A — > 1, which fact 
alone favors choosing A close to unity. One manifesta- 
tion of discretization error is a modification of the effec- 
tive coupling between impurity and bosonic bath [con- 
duction band] . Essentially exact results, if desired, can 
be obtained by performing a A — > 1 + extrapolation from 
A > 1 results. Fortunately, critical exponents are often 
found to be extremely insensitive to the chosen A. Larger 
A (e.g., A = 9) calculations, which retain fewer states and 
thus require less computational effort, can then be used 
to extract the critical properties of the model. However, 
the weak A-dependence of critical exponents cannot of 
course be assumed at the outset. 

Finally, Eqs. 1)3. 8fl and i)3. 9fl are mapped exactly onto 
tight-binding Hamiltonians via a discretized Lanczos 
procedure. 50 Defining 



E 



2 _ 
u nm UJ 'm — e m 



(3.12) 



fna — y ] ( v nrn c rncr + v nm c m(j)i (3.13) 



where Uo m = B m /B and v^ rn — l/\/2, such that 
[°n, &„/] = S n:7l ' and {f na , fl^ a ,} = &<j,<j' , we obtain 

H B = u Q [eX&n + tn(&J,&n-i + , (3-14) 



with oj m given by Eq. I|3.10|l . For a symmetric fermionic 
band density of states p(e) = p(—e), e„ = for all n. For 
the specific case Eq. (|2.7|l considered in the following, we 
use the closed-form algebraic expression for r„ as derived 
by Wilson»ii The transformed Hamiltonian to be solved, 
H = H F + H B + H int [Eqs. (2.9), JXTHl and (TTTBTl l. 
describes two semi-infinite chains each coupled at the end 
site only to a spin-i magnetic moment. A schematic of 
the chain form of the Hamiltonian is depicted in Fig. 0J. 



B. Iterative solution 

The on-site energies e„ and hopping coefficients r„ for 
the fermionic tight-binding chain decrease like A"™/ 2 for 
large chain length n. The hopping coefficients for the 
bosonic chain, e n and t n , drop off faster — as A - ' 1 for 
large n — arising because the spectral density B(lj) is 
nonzero for lu > only. The latter fact has been dis- 
cussed in the context of a pure bosonic NRG approach 
to the spin-boson model^^ 5 " and there allows a conver- 
gent approximation using a single bosonic chain of fi- 
nite length. While the BFKM mapping is essentially a 
combination of the pure fermionic and pure bosonic ap- 
proaches, the iterative procedure described below is an 
important novel feature of our approach, and one that 
does not a priori guarantee a satisfactory description of 
the quantum-critical physics of the model. 

In employing the NRG for the BFKM we must ensure 
that at each step we treat only fermions and bosons of 
the same energy scale. As for the pure fermionic case, we 
take D as the energy unit and write the full Hamiltonian 
as the limit Eq. I|3.1|l . where here N refers to the highest 
labelled fermionic site or iteration number. Crucially, the 
following recursion relation, relating to i?jv-i(M — 
1), adds site M to the bosonic chain only if its tight- 
binding coefficients are comparable to the energy scale 
{Ar N / 2 ) set by the fermionic site N added at the current 
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iteration, N. Thus, 

a 

+ Q- 1 A JV / 2 ^r J v(/L/jv-i, CT + H.c.) 

(7 

+ Qoa-W'QteM - ck- N ' 2 ) \e M tf M b M 



JVtr 



tjVf(&M&M-l + H.C.) 



(3.17) 



where u)o = ^q/D, £m = naax(ejw;ijvf) and c is a 
threshold parameter that we typically assign the value 
c = A~ 1//4 . The Hamiltonian for iteration N = 0, con- 
taining the couplings between the impurity and the end 
site of each chain, is given by 

aH = ^2 poJoS ■ /o er <w/oa-' + X! e o/o<r/o<T 

crcr' cr 

+ u ^p^-S z {bl + b ) +u e blb , (3.18) 

with natural dimensionless parameters PqJq, Bq = 
(Kq9o) 2 and ojq = ujq/D to be varied for a given exponent 
s in the bath spectral function B(lo) [Eq. (|2.8|) ]. 

Since both e n and i„ decay as A - " for large chain 
length n, the iterative procedure Eq. I|3.17|l rapidly 
amounts to extending the bosonic chain at every second 
iteration (i.e., for every second site added to the fermionic 
chain). Thus without any change in the essential phys- 
ical properties, it is possible, and proves convenient, to 
modify Eq. (|3.17|) to force the extension of the bosonic 
chain to occur only at a given iteration parity. 

We also note in passing that an alternative NRG 
scheme for the BFKM might extend both chains at each 
iteration but employ a different discretization parameter 
Kb for bosonic energies related to that for fermionic en- 
ergies Ap by Ap = A 2 B (Rei. IFli) . Though such a scheme 
has not yet been implemented, we believe it does not con- 
fer any significant advantages over the approach taken in 
this work. 

The eigenstates {\E)n-i} of H^^i are used to con- 
struct the basis states for iteration N: 



\l,E) N = \E) N _ 1 ®\l) N . 



(3.19) 



where {|Z)jv} are the states of the added site(s). If the 
fermionic chain alone is extended then |Z)jv = |i)jV) onc 
of the four possible states of the added fermionic site. 
If both a fermionic and a bosonic site are added, then 
\1)n = \i)N ® Wm) where{|nA./)} is an eigenstate of 
b\fbM- In contrast to fermionic sites, the number of 
bosons per site of the bosonic chain is in principle un- 
limited. In practice, the bosonic Hilbert space must be 
truncated, restricting the maximum number of bosons 
per site of the bosonic chain to a finite number Nb, typi- 
cally Ni, < 16 and thus including only the lowest Nb + 1 
eigenstates of b' M bfri. We have systematically studied 
the effects of this truncation on various properties of the 



BFKM. In particular, we find that universal critical prop- 
erties rapidly converge with increasing Nb — typically 
8 < Nb < 16 suffices — providing a very satisfactory 
description of the QCP for the model. This is discussed 
further in Sec. IIVI 

It follows from Eqs. (|3.17|) and (|3.19|l that the matrix 
elements of Hn, n(V , E'\Hn\1, E)n, can be evaluated 
in terms of the eigenstates and eigenenergies of Hm-i 
and the matrix elements N{i'\f nV) n and, if appropriate, 
m (n' Ib^jln) m ■ The numerical effort required to diago- 
nalize Hn is considerably reduced by taking advantage 
of various symmetries of the Hamiltonian, as discussed, 
e.g., in Refs. 1491 and l34l With Hn diagonalized, the ma- 
trix elements of H^+i are formed in the basis Eq. (|3.19() 
with TV — » TV + 1 and the procedure is repeated to obtain 
the many-body spectrum on progressively lower energy 
scales. 

The need for an iterative approach arises because (even 
after bosonic truncation) the size of the basis rapidly be- 
comes unmanageable. For the regular Kondo model it 
grows by a factor of 4 at each iteration, reflecting the 
possible occupancies of the added fermionic site. Even in 
that relatively simple application, finite computational 
power demands that the basis be truncated after only 
a few iterations, saving only the lowest N s many-body 
eigenstates of H^-i to be used to construct the basis 
in which to diagonalize -ff/v. The problem is of course 
even more acute in the BFKM: for each iteration that a 
bosonic site is added to the chain the basis grows by a 
factor of4(iV 6 + l). 

By retaining only the lowest N s eigenstates from one 
iteration to the next, the NRG method gives a sequence 
of truncated Hamiltonians from which thermodynamic 
[dynamical] properties may be calculated. For each iter- 
ation N, the error due to truncation is minimized in a 
temperature [frequency] window around a characteristic 
energy scale 



T N = aA- N / 2 /$, 



(3.20) 



where (3 < 1 is a dimensionless parameter— that gener- 
ally decreases with increasing N s . Properties at temper- 
ature [frequency] scales smaller than Tjy will be more ac- 
curately calculated at subsequent iterations, while those 
at higher scales have been more accurately calculated at 
a previous iteration. See, e.g., Refs. 1491 andl5^ for further 
discussion. 

In the following we turn to results arising from the 
above approach, beginning with the case of sub-Ohmic 
dissipation. 



IV. RESULTS: SUB-OHMIC DISSIPATION 

We now consider results for the Bose-Fermi Kondo 
model with sub-Ohmic dissipation, i.e., s < 1 in Eq. 
I|2.14[) . The Ohmic case s = 1 is discussed separately 
in Sec. For < s < 1 the model exhibits a continu- 
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FIG. 2: Schematic of the NRG flow of effective couplings for 
(a) < s < 1 and (b) — 1 < s < 0. Trajectories describe 
the flow with increasing N (or decreasing Tjv) of the effective 
J± and g that enter Hn- For < a < 1 an unstable critical 
fixed point (solid circle) lies on the separatrix (dashed line) 
that specifies the phase boundary between the stable Kondo 
(K) and localized (L) phases in the J±-g plane. The Kondo 
fixed point corresponds to an effective J± = J z = oo and 
g = 0, while the localized fixed point has J± = 0, J z finite, 
and g = oo. For — f < s < the localized fixed point is 
reached for all g =fc 0. See text for further discussion. 



ous impurity quantum phase transition between Kondo- 
screened and unscreened-local-moment phases. After dis- 
cussing the NRG flows, fixed points and phase bound- 
aries, we demonstrate that a low-energy scale T„, char- 
acterizing the crossover from critical to stable fixed 
point behavior, vanishes in a power-law fashion at the 
quantum-critical point. Local magnetic properties con- 
sidered in Sec. IIVDI show distinct behaviors in either 
phase and at the QCP itself, from which various criti- 
cal exponents can be extracted. The interacting nature 
of the QCP is evident from hyperscaling of the criticial 
exponents, including o;/T-scaling behavior in the dynam- 
ics. We conclude the discussion of the sub-Ohmic case 
by considering the local spectrum of single-particle exci- 
tations for the related Bose-Fermi Anderson model. The 
destruction of the central Kondo resonance at the QCP 
is clearly observed. 

For notational compactness we write henceforth 
p J = J and K g = g. We also set the bosonic bath 
cut-off u>o = 1 in our calculations^ although we em- 
phasize that our approach is in no way limited to this 
particular choice. All energies, frequencies, and temper- 
atures are expressed as multiples of D = 1. For a given 
bath exponent s, we typically work at fixed J and vary g 
to tune the system through a quantum phase transition. 
However, in calculations of the single-particle spectrum, 
we specify Bq = g 2 . Critical properties are found to be 
insensitive to the NRG discretization parameter A and 
adequately converged for bosonic truncation parameter 
8 < A;, < 12. Unless stated otherwise, results presented 
were obtained using A = 9, Nb = 8 and N s = 500. 



A. NRG flow and fixed points 

Figure Ufa) shows a schematic of the NRG flow of ef- 
fective couplings found for < s < 1. There are two 
stable fixed points of the flow: the familiar Kondo fixed 
point (K), where the impurity spin is screened by conduc- 
tion electrons and the single-particle spectrum contains a 
Kondo resonance (see Sec. lIVEl ): and a "localized" fixed 
point (L), where the impurity dynamics are controlled by 
coupling to the dissipative bath. A third fixed point — 
the unstable, critical fixed point (where Kondo screen- 
ing is critical) — lies on the separatrix g — g c {J) which 
specifies the phase boundary in the J-g plane (see Sec. 

hvbv 

As an example, Fig.|21shows the NRG flow of the low- 
lying eigenenergies of Hn as a function of even iteration 
number N, for fixed Kondo coupling J and five different 
couplings g to a sub-Ohmic bosonic bath with s = 0.2. 
Increasing N corresponds to considering the problem on 
a sequence of exponentially decreasing energy scales [Eq. 
(|3.20|) ]. The bosonic chain was extended at each even 
iteration starting at N = 4. Figure EJa) shows the flow 
for g — 0, where the impurity spin and bosonic bath 
are decoupled. The fixed-point spectrum reached, with 
effective J infinite, corresponds to that of the regular 
Kondo modeli^ plus a chain of free bosons. The resulting 
fixed-point Hamiltonian, H* K = R 2 [HjA, thus consists 
of a chain of free bosons and a chain of free fermions 
described by 

JV 

H { F L) N = a^A^ 2 tnfUn,* + r n (flf n - ha + H.c). 

n—L 

(4.1) 

with L = 1, reflecting the fact that the local moment 
and conduction electron state at the impurity (/o) are so 
strongly coupled — via an infinite Kondo interaction J 
- that those degrees of freedom are frozen out. Each 



TABLE I: Low-lying eigenstates of the even-N fixed-point 
Hamiltonian obtained for s — 0.2. The {77* } denote the 
limiting eigenvalues of Hp N [Eq. 14. H i corresponding to the 
regular Kondo model fixed-point spectrum, while {£* } are the 
limiting eigenvalues for a chain of free bosons. 
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FIG. 3: Low-lying energy levels of Hn as a function of even iteration number N, for bosonic bath exponent s = 0.2, Kondo 
coupling J = 0.5 and the values of g shown, with g c — 0.567. The bosonic chain was extended for even N only, as discussed in 
Sec. 1111 Bl (a) Zero coupling to the bosonic bath, (b) For g close to g c in the Kondo phase, the quantum-critical fixed-point 
regime [shown in (c)] is accessed for intermediate N (here 10 < N < 20) before the levels flow to the Kondo-phase fixed-point 
structure (for N > 40). (d) and (e) In the localized phase, the levels cross over from the critical fixed-point [10 < N < 20 in 
(d)] to the localized fixed-point structure observed for large N . The crossover in the flows defines a low-energy scale T* that 
vanishes as g — > g c . 




N (even iterations) 



FIG. 4: Interpretation of the low-lying energy levels of the 
localized fixed-point Hamiltonian H^. The left-hand portion 
of the figure shows Hn as a function of even iteration number 
N, for s = 0.2, J = 0.5, and g — 0.7. The spectrum can 
be decomposed into the product of pure-fermionic and pure- 
bosonic spectra. Plotted on the right-hand side are the levels 
listed in Table HJ with pure-fermionic excitations shown as 
dotted lines, pure-bosonic excitations as thick dashed lines, 
and mixed excitations shown as thin dashed lines. 



chain problem can be numerically diagonalized to yield 
a set of single-particle levels that does not depend on 
N for N ^> 1 (apart from an odd-even alternation for 
the fermionic case). A decomposition of the lowest lying 
eigenstates of the even-A^ Kondo fixed point is provided 
in Table HI 



For all g < g c , the effective bosonic coupling flows to 
zero [see Fig.|2Ia)], while the effective J flows to infinity, 
and the fixed-point level structure reached for large itera- 
tion number N is identical to that for g = 0. This is seen 
clearly in Fig.[3Jb). However, for g close to g c , the flow 
first approaches an unstable or critical fixed point. This 
is also clearly observed in Fig.^b) for 10 < N < 20, be- 
fore the flow crosses over to the stable Kondo fixed point. 
The low-energy scale set by the crossover from unstable 
to stable fixed point vanishes as g — > g c and is discussed 
further in Sec. IIV CI The level structure for the critical 
fixed point is shown in Fig. Etc), obtained for g ~ g c . 

In the localized phase g > g Cl the spectrum of the fixed- 
point Hamiltonian iJ£ can also be decomposed into a 
product of pure-bosonic and pure-fermionic spectra. Fig- 
ureHJshows NRG level flows for g = 0.7 and otherwise the 
same parameters as Fig. [21 Localized-fixed-point behav- 
ior is there seen for N > 15, and the fixed-point spectrum 
can be decomposed into pure-fermionic excitations, pure- 
bosonic excitations, and appropriate linear combinations 
of the two, as shown in Table ITT1 The bosonic spectrum 
is identical to that of the localized fixed point of the 
corresponding sub-Ohmic spin-boson model, while the 
fermionic spectrum corresponds to that of an anisotropic 
Kondo model in which J± has renormalized to zero [see 
Fig. Ufa)] . The fixed-point value of J z is finite for g > g c 
and diverges as g — » gj according to 

J* z K(g-g c )-P, (4.2) 

where /3 is the magnetic critical exponent defined in Sec. 
IIVDI This behavior is demonstrated in Fig. [S] which 
shows J* versus g — g c for J = 0.5 and two values of 
the bosonic bath exponent. In both cases the values of 
the exponent (3 agree with those obtained in Sec. lIVlJl to 
within 2%. 

As s — > + the critical fixed point merges with the 
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FIG. 5: Effective coupling J* at the localized fixed point for 
the bath exponents specified in the legend and J = 0.5. J% 
exhibits a power-law divergence as <?—><? t [see Eq. 14.21 . 
The values of the exponent (5 agree with those obtained in 
Sec. lTVDl to within 2%. 
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FIG. 6: Critical coupling g%, separating Kondo and localized 
phases, versus bosonic bath exponent s for fixed Kondo in- 
teraction J = 0.5. The inset shows the Tk dependence of 



dashed lines are fits to the form gl cx T K 
for Tk <C 1. See text for discussion. 



which is followed 



Kondo fixed point. For any — 1 < s < 0, the critical 
coupling g c {s) = 0, i.e., the system is localized for any 
nonzero g. For small g the flow approaches the Kondo 
fixed point discussed above [shown schematically in Fig. 
|2Jb)], before crossing over to the localized fixed point on 
an energy scale T* that vanishes as g — > 0. For s = 1 



TABLE II: Low-lying eigenstates of the even-iV fixed-point 
Hamiltonian iJ£, which can be decomoposed into the prod- 
uct of a pure bosonic and pure fermionic spectrum. The 
former, {w*}, is that of the localized fixed point for the 
spin-boson model and is independent of g > g c . The lat- 
ter {r/*} correspond to a pure- fermionic anisotropic Kondo 
model: J*(g = 0.7) ~ 2.772 and = 0. The fixed-point 
value J* diverges as g — * 
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the QPT is Kosterlitz-Thouless-like, considered in Sec. 
IVl while for s > 1 the Kondo fixed point is reached re- 
gardless of the strength of coupling to the bath. 



B. Critical coupling 

The critical coupling g c , which marks the phase bound- 
ary between Kondo-screened and localized phases for 
< s < 1 , is shown in Fig. as a function of bath expo- 
nent s. For s — > + , g 2 c (s) vanishes linearly in s. The line 
of continuous quantum phase transitions for < s < 1 
terminates at the Kosterlitz-Thouless transition point at 
s = 1. 

Figure [7| demonstrates, for representative bath expo- 
nents s = 0.2 and 0.8, that the critical couplings con- 
verge rapidly upon increasing the truncation parameters 
N s and 2V&, denoting respectively the number of states re- 
tained from one NRG iteration to the next and the max- 
imum number of bosons allowed per site of the bosonic 
chain. In percentage terms, g c converges with increasing 
Nb more rapidly for smaller s, but converges with increas- 
ing N s more rapidly for larger s. However, in all cases 
the choices N s ~ 500 and 8 < Nf, < 12 are sufficient for 
the determination of phase boundaries for A = 9. In Sec. 
HVDI we show that critical exponents are also adequately 
converged for these parameters, a result of crucial im- 
portance in judging the success of the approach, and one 
that is not a priori guaranteed. 

Of course, deep inside the localized phase it is to be 
expected that an approach using the lowest Nb + 1 eigen- 
states of b^jbM as basis states for each bosonic site M, 
will ultimately fail due to a divergence of the mean site- 
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0.578 




FIG. 7: Dependence of critical coupling g c on NRG truncation 
parameters JVj, and N s for J = 0.5, NRG discretization A = 9, 
and bath exponents s = 0.2 and s = 0.8. The first column 
demonstrates rapid convergence in Nb for fixed N s — 500. 
The second column shows convergence in N s for fixed N b — 8. 

occupancy. In a recent study of the SBM using a bosonic 
NRG approach, 44-45 Bulla et al. have employed a bosonic 
basis of displaced oscillator states implemented within 
an alternative "star" formulation of the NRG. While this 
approach describes well the localized fixed point, it fails 
to capture correctly both the flow to the delocalized fixed 
point and the critical properties of primary interest 4 s - 

The preceding comments are well-illustrated by Fig. 
IH1 which considers the operator for the total number of 
bosons at iteration N, 

M(N) 

B N = b\h, (4.3) 

i 

and its expectation value evaluated at the characteristic 
temperature T N , Eq. IpOUJ) . Here, M(N) denotes the 
highest labelled bosonic site at iteration 7Y. Figure |S| 
shows (-B20) versus bosonic truncation parameter Nf, for 
both Kondo (g < g c ) and localized (g > g c ) phases and at 
the QCP (g = g c ) itself. In the Kondo phase and at the 
QCP, expectation values converge rapidly in Nb, indicat- 
ing that the lowest eigenstates of Hm are well-described 
using only the lowest Nb + 1 eigenstates of b' M bM as ba- 
sis states for each bosonic site M. It is clear that no 
such convergence is achieved inside the localized phase. 
The difference between the saturated values of (-B20) f° r 
critical and delocalized phases is smaller for lower s, con- 
sistent with the merging of the critical and delocalized 
fixed points as s — > (see Fig. [2}. 

The aim of the present work is to provide a thorough 
numerical account of the critical behavior of the BFKM. 
As is demonstrated explicitly in Sec. IIV Dl the approach 



FIG. 8: Dependence of (-B20), defined in text and evaluated 
at the characteristic temperature scale of iteration N = 20, 
on bosonic truncation parameter Nb for (a) s — 0.2 and (b) 
s — 0.8. Data were obtained for J = 0.5 and the bosonic 
couplings g shown in the legend. Convergence with respect 
to Nb is reached rapidly for g < g c . 

described in Sec. IIIII provides such a description. We em- 
phasize that critical properties are well-captured regard- 
less of the phase — Kondo or localized — from which 
the critical point is approached. 

We now return to a discussion of the phase boundary 
between Kondo and localized regimes of the model. Re- 
ducing J, and hence the Kondo scale Tk of the pure- 
fermionic (g = 0) problem, squeezes out the Kondo- 
screened phase of the BFKM; specifically, we find 

g 2 c {s) = c(s) T]f s for T K « 1, (4.4) 

with c(s) an s-dependent constant that vanishes both as 
s — > + and s — > 1~, in agreement with perturbative (in 
e = 1 — s) RG calculations iii 

The form of Eq. I|4.4|l is demonstrated in the in- 
set to Fig. |B| which shows g\ versus Tk on a loga- 
rithmic scale for s = 0.2 and 0.8. Here, we define 
Tk = Dexp[—A(A)/J]. The A-dependent constant A(A) 
tends to correct the effects of discretization, which re- 
duce the effective coupling between impurity and band, 
and leads to the most faithful description of the impu- 
rity problem^ The slopes of the linear fits, shown as 
dashed lines, reveal that the exponent is 1 — s to high 
accuracy, which relation holds for all < s < 1. We 
note that this exponent follows from simple dimensional 
arguments. When all other energy scales in the prob- 
lem are much smaller than the cut-offs, the problem may 
be rescaled in terms of Tk, with a single dimcnsionlcss 
parameter g 2 (TK /wo) 1_s '. 

If the Kondo scale Tk is recognized as a (renormal- 
ized) tunneling amplitude between the two spin states of 
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9: (a) Power-law behavior of the crossover scale T, 



N */ 2 in the vicinity of the QCP. Data are shown for s = 0.5, 
= 0.5 (solid symbols) and J = 0.1 (open symbols), for both 



the Kondo (circles) and localized (squares) phases with NRG 
discretization A = 9. Kondo-phase data obtained for s = 0.5, 
J = 0.5, and A = 3 are shown as stars. The correlation- 
length exponent v(s) is universal and is independent of J, of 
the phase from which the QCP is approached, and of A. (b) 
Schematic crossover diagram valid near the QCP for < s < 
1. The dashed line shows T*, the scale for crossover between 
critical and stable fixed-point behavior in static or dynamical 
quantities. The horizontal dashed line marks Tk, the Kondo 
scale of the pure-fermionic problem, which serves as a high- 
energy cutoff for the quantum-critical behavior. 
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FIG. 10: Dependence of the correlation length exponent v(s) 
on bath exponent s. Results for the spin boson model are 
plotted as crosses. For s — » 0, v^ 1 ~ s (fit shown as dotted 
line) while v^ 1 vanishes as — s) for s-»l (fit shown as 

dashed line). 



the impurity, then Eq. (|4.4H follows from the sub-Ohmic 
spin-boson model result f*4 where it is found that the crit- 
ical dissipation strength a c oc A 1_s for A <C luq. Such 
an interpretation follows, of course, from the bosoniza- 
tion discussed in Sec. Ill Bl The sub-Ohmic BFKM can 
be mapped to a spin-boson model with a sub-Ohmic bath 
plus an Ohmic contribution arising from bosonization 
of the Kondo part of the Hamiltonian. While the sub- 
Ohmic behavior dominates the bath spectral function at 
low energies, and hence determines the universal critical 
properties, the Ohmic contribution acts to renormalizc 
the bare couplings, as discussed more generally in Ref. 
Specifically J — > Tk, leading to an effective SBM with 
A oc Tk and a sub-Ohmic bath. 

C. Crossover scale 

We now extract the correlation length exponent v for 
the BFKM, considering bath exponents in the range 
< s < 1. It is a general feature of a continuous quantum 
phase transition^ that a characteristic energy scale for 
fluctuations above the ground state vanishes as g — > g c 
with a universal exponent vz. The exponent v describes 
the divergence of the spatial correlation length at the crit- 
ical coupling, and the dynamical exponent z specifies the 
effective spatial dimensionality of each time dimension. 
Note that z = 1 for impurity models. 

Close to the critical coupling g c , the crossover from the 
unstable quantum-critical regime to one or other of the 
stable regimes, defines an energy scale T* that vanishes 
precisely at g c in a fashion governed by the correlation- 



length exponent: 

T m x\g-g c \ v . (4.5) 

T^T 1 sets the time-scale for the decay of the local- moment 
fluctuations. Close to the QCP the characteristic energy 
scale in the problem is , while at the QCP itself the only 
relevant energy scale is the temperature, as discussed in 
Sec. lIVD2l in the context of uj/T scaling. The crossover 
scale T* may be determined directly from the NRG many- 
body spectrum, via T* oc A' N */ 2 . Here, N* is the NRG 
iteration number at which a selected energy level crosses 
over from the critical fixed point to a stable fixed point. 
For example, for the lowest nonzero energy level in Fig. 
EJb), the clear crossover from critical (E c ~ 0.246) to 
Kondo regime (E K — 0.322) occurs for N* ~ 27.2, with 
here defined such that En^ — \ {E C + Ek)- Alterna- 
tively, T* may be extracted from crossovers seen directly 
in properties such as the local static or dynamical suscep- 
tibility or the single-particle spectrum. Good agreement 
is found in all cases to within the numerical uncertainty, 
which is typically less than 1%. 

Fig. Efa) shows T* oc iV* versus \g — g c \ for J = 0.5, 
s = 0.5 and a range of couplings g approaching g c ~ 1.22 
from either phase. The form Eq. (|4.5() holds well over 
several orders of magnitude, as reflected in the linear be- 
havior of the data on a logarithmic scale, and any devi- 
ations for small \g c — g\ are removed by a more accurate 
determination of g c . The correlation-length exponent v 
is proportional to the slope of the lines and it is there- 
fore clear to see that v is independent of Kondo coupling 
J and the phase (Kondo or localized) from which the 



13 



QCP is approached. The exponent is universal, as are all 
others discussed below, and depends only on the bosonic 
bath exponent s, which thereby plays the role of a dimen- 
sion. The insensitivity to A is also demonstrated in Fig. 
IHJa), justifying our choice to work with A = 9 (and hence 
fewer retained states) for extracting critical exponents. 

The behavior shown in Fig. Efa) is representative of 
that for all s in the range < s < 1. A generic crossover 
diagram is shown in Fig. EJb), valid in the vicinity of 
the QCP and for T (and u>) below a nonuniversal cutoff 
scale (above which free-local- moment behavior obtains). 
For T < T*, behavior characteristic of either the Kondo 
(g < g c ) or the localized {g > g c ) phase is observed, while 
for T > T* behavior characteristic of the QCP itself ob- 
tains. At the QCP, where T* = 0, the latter behavior per- 
sists right down to T = 0. This results in the familiar pic- 
ture for a quantum phase transition of quantum-critical 
behavior emanating from a singular point at T — and 
influencing a large part of the phase diagram in the g- 
T (or g-u) plane. Such a picture has emerged in recent 
years in a number of quantum impurity problems^ 

The dependence of v on the continuous parameter (di- 
mension) s is shown in Fig. llOf ah For both s — > + and 
s — > 1 _ the correlation length exponent diverges, which 
reflects the qualitative changes occuring at these critical 
dimensions, as described above. For small s, 



v{s) ~ 1/s, 



(4.6) 



which result holds reasonably well for s = 0.4, and be- 
comes asymptotically exact as s — > + [see dotted line fit 
in Fig. flUt a)]- For s — > 1 the exponent approaches 



v{a) ~ 1/ V2(l - a) + C, 



(4.7) 



with C ~ 0.4, which fit is shown as the dashed line in Fig. 
HOf a). Most importantly, we demonstrate that the corre- 
lation length exponents of this Ising anisotropic BFKM 
are in essentially exact agreement with those of the spin- 
boson model obtained in Refs . 1441 and l48l This is true for 
all exponents calculated, providing numerical confirma- 
tion that the quantum-critical points of the two models 
belong to the same universality class — an equivalence 
implied by the bosonization treatment discussed in Sec. 

una 

The results v(s) ~ 1/s and gl/T-^ oc s as s — > [Eqs. 
(|4.6() and (|4.4() respectively] imply a form for the vanish- 
ing of the crossover scale T* at s — 0. Here the system 
is always localized for any g > and we predict that T„, 
describing the direct crossover from Kondo to localized 
behavior, vanishes approaching g = g c — as 



lnT* oc 



g 2 



for s = 0. 



(4.8) 



This is confirmed numerically in Fig. 1111 which shows 
data for J = 0.5 rescaled according to Eq. I|4.8|l . The 
dashed line shows a linear fit for g 2 <§; 1. The inset 
demonstrates that the slope of the fit is proportional to 




FIG. 11: Crossover scale for s = vanishing as Eq. 1)4. 8|l 
for g — > 0. The data in the main figure were obtained for 
J = 0.5. A linear fit to the data points at lowest g 2 is shown 
as a dashed line. The inset demonstrates that the slope m of 
the fit varies as 1/Tk — see text. 



T" 1 , where T K = Dexp[-A(A)/J]. For the sub-Ohmic 
spin-boson model, identifying A oc and a oc g 2 leads 
to lnT* oc (—A/a) for s — 0, which result has also been 
obtained via a perturbative RG analysis. 48 



D. Local magnetic response 

The local susceptibility, xiocM = Xi oc M + *Xk,c (<*>)> 
is defined as 



XiocM = i 



f ([S z (*),^(0)]) (4.9) 



with XiocC^ = 0) being the static response 



Xloc(0) 



d(S z 



dh 



h=0 



lim 

h^O h 



(4.10) 



to a local magnetic field h acting only at the impurity site 
through an additional term in the Hamiltonian AH = 
hS z (setting gfiB = !)• 

In the present section we focus on the sub-Ohmic case 
with < s < 1. The case of Ohmic dissipation, s = 1, is 
discussed separately in Sec.E] In Sec. UVDD we consider 
the local static susceptibility \\oc(T; u> = 0) calculated by 
evaluating the right-hand side of Eq. I|4.10|l directly for 
small fields in the range 10 -12 < h < 10~ 8 . We ex- 
tract various critical exponents for the transition which 
are shown to satisfy hyperscaling relations characteristic 
of an interacting quantum-critical point. We calculate 
the dynamical local susceptibility in Sec. IIV D 21 which 
obeys lo/T scaling in the vicinity of the QCP. We reiter- 
ate that all results presented are converged with respect 
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FIG. 12: Static local susceptibility xioc(T',(^ = 0) versus T 
for s = 0.5, poJ = 0.5 and different g, where g c — 1.21583. 
Xloc(0,0) = for g < g c . At the critical coupling we find 
X.\oc(T;0) cx T~ x with x = 0.500(2). The inset shows the 
vanishing of Txi oc (0:0) ^ or ^ fl^- Analogous data for s 
0.2 are shown in Fig. 2 of Ref. f 
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FIG. 13: Dependence of the critical exponent x, as defined in 
Eq. 14. 161 . on bath exponent s in the range < s < 1. The 
data were obtained for Kondo coupling J — 0.5, discretization 
parameter A = 9, and (a) N s = 500, (b) N b = 12. The dotted 
lines show x — s. 



to the truncation parameters Nf, and N s . For the data 
discussed, Nb — 8 and N s = 500 are sufficient to converge 
critical exponents to the required accuracy — see, e.g., 
Fig. 

1. Static susceptibility 

The Kondo and localized phases are readily distin- 
guished by the limiting low-temperature behavior of 
X\oc{T\lo = 0). If a local moment M\ oc = (S z ) exists 
at the impurity site for T — > 0, then 

Xlo c(T; u = 0) = for T « T*. (4.11) 

We find such behavior for g > g c , as shown in Fig. ITU 

For g 3> g c the residual moment Mi oc (0) approaches 
the value for a completely decoupled spin, Mi oc (0) 2 = 
S(S + l)/3, as seen clearly in the inset to Fig.^1 For < 
s < 1, the low-temperature limit of Txioc(.9 > 9c,T;uj — 
0) vanishes continuously at the QCP, viz 

lim T ^ rxioc(g > 9c, T;u = 0)cx{g- g c )\ (4.12) 

as is also clear from the inset to Fig.^21 Equations (|4.11l) 
and H4.12fl imply that 

M loc ( 5 >« ?c ;T = 0)cx( 3 -. 9c )' 3 (4.13) 

with (3 = A/2. 

For g < <7 C , by contrast, we find the characteristic sig- 
nature of a screened ground state with quenched impurity 



degrees of freedom (as exists ubiquitously in the regular 
Kondo model): 

Xioc(T;w = 0) = const, for T<T», (4.14) 

which behavior is also seen clearly in Fig. For < 
s < 1 it is found that 

X!oc(3 < 9c\ T = 0, u; = 0) <x (g c ~ gY\ (4.15) 

signalling a collapse of the Kondo effect at g = g c (see 
Fig. I14|) discussed further below. 

At the QCP itself, g = g c , we find the anomalous be- 
havior 

Xioc{T;uJ = 0,g = g c )^T- x . (4.16) 
To high accuracy we find 

x = s (4.17) 

for < s < 1, confirming the prediction of e-expansion 
studiesiii^ that x = s to all orders in e = 1 — s. The re- 
sult Eq. (|4.17|l has also recently been reported for the 
sub-Ohmic spin-boson model4^& We note in passing 
that since the total spin quantum number is not con- 
served in the BFKM, the impurity contribution to the 
total susceptibility Ximp also acquires an anomalous ex- 
ponent at criticality. Again, we find XimpC? 1 ) oc T~ s . 

For T <C Tk the local static susceptibility obeys the 
scaling form 



TkX\oc{T\u = 0) = 




(4.18) 
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FIG. 14: Order parameter Mi oc (g > g c ;T = 0) and renormal- 
ized Kondo temperature Xloc(<7 < 9c\T = 0,u = 0) _1 versus 
bosonic coupling g, for s = 0.2 and J = 0.5. Both quantities 
vanish continuously for g — > g c . The inset shows fits (dashed 
lines) for small \g — g c \. See text for discussion. 



FIG. 15: Dependence of residual moment M\ oc {g — g c ,T = 0) 
on local magnetic field h, for J = 0.5. The dashed lines are 
fits to the data, as discussed in the text. The inset shows the 
data on an absolute scale. 



where 4> StP is a universal function andp = K or L denotes 
the phase. For z» 1, <j) s ^ p (z) oc z~ s and TkX\oc{T) oc 
{T/T K )~ S is independent of T*. 

Before preceding we illustrate the convergence of the 
critical exponents in the NRG truncation parameters Nf, 
and N s , which denote respectively the maximum number 
of bosons allowed per site of the bosonic chain and the 
maximum number of states retained from one NRG iter- 
ation to the next. The behavior of the critical exponent x 
is typical and is shown in Fig. 1131 as a function of bosonic 
bath exponent s. Figure H~3T a) shows data obtained for 
Nf, = 4, 8 and 12, each converged with respect to N s , the 
number of states retained. Upon increasing Nb, the ex- 
ponent x is found to approach x — s [dashed line in Fig. 
I13f a)] for all s, the convergence being faster for small 
s. For this and all other critical exponents calculated 
we find that by 8 < Nb < 12 sufficient convergence is 
achieved. Figure H^T b'l illustrates the dependence of x on 
the number of states retained N s , calculated for Nb = 12 
and three values of s. The exponent converges rapidly 
with increasing N s , approaching x — s (dashed lines) in 
each case. In summary, we find that by Nb — 12 and 
N s = 500 the critical exponent x is converged to within 
1% of x = s for all s examined in the range < s < 1. 
The observation that fewer bosons are required to achieve 
convergence for s <C 1 is consistent with the merging of 
the critical and delocalized fixed points as s — + 0. 

We now return to a discussion of the various critical ex- 
ponents introduced above, for which similar convergence 
is observed. The vanishing as g — > of the residual mo- 
ment M\ oc (T = 0), which serves as the order parameter 
for the problem, is shown in Fig. 1141 (squares) . This pic- 



ture remains qualitatively valid for all bath exponents in 
the range < s < 1. For the specific case s = 0.2 shown, 
the exponent defined in Eq. (|4~T3|) is /3 = \\ ~ 1.99(2), 
with the estimated numerical uncertainty in the last digit 
given in parentheses. The fit is shown as a dashed line in 
the inset to Fig. El Also shown in Fig. El is the vanish- 
ing of x\oc{g < g c ;T = 0, w = 0) for g —> g~ described by 
exponent 7 [Eq. 14.15) . xioc(0, O)^ 1 measures the Kondo 
temperature, which is strongly renormalized due to the 
competiting coupling to the dissipative bosonic bath, and 
its vanishing as g — > g~ signals criticality in the Kondo 
screening. For s = 0.2, 7 ~ 1.00(1), and the dotted line 
in the inset shows a linear fit to the g — ► g~ data points. 

The dependence of M\ oc (h; g = g c , T = 0) on the local 
field h defines the magnetic critical exponent 5: 

M loc (h;g = g c ,T = 0) oc \h\? . (4.19) 

FigurelTSlshows data obtained for J = 0.5 and g = g c for 
both s = 0.2 and s — 0.8; the inset shows the data on an 
absolute scale, while the main figure reveals the power- 
law behavior that is observed over several orders of mag- 
nitude. The dashed lines are fits to the h < 10 -3 data, 
of the form of Eq. (O^ . with 1/6 ~ 0.665(5) [0.106(2)] 
for s = 0.2 [0.8]. 

Hyperscaling relations among the various critical expo- 
nents may be deduced starting from the following scaling 
ansatz for the critical part of the free energy, which as- 
sumes that the fixed point is interacting: 

F CYit = Tf(\g-g c \/Ti,\h\/T b ). (4.20) 
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FIG. 16: Behavior of critical exponents 0, 7 and 1/5 as a func- 
tion of bosonic bath exponent s in the range < s < 1. The 
symbols correspond to directly calculated exponents, while 
the solid lines were obtained from hyperscaling relations [Eqs. 
(4.23]) using numerical results for x(s) and v(s). Hyperscal- 
ing relations are obeyed to within numerical accuracy for all 
s considered in the range < s < 1. 



Using the standard thermodynamic relations 

<9-Fcrit 



loc 



dh ' 



X loc(l> = 0) = — 



(4.21) 
(4.22) 



dh 2 ' 

it readily follows that there are only two independent 
exponents, i.e., 



1 + x 

d = 1 ' 

1 — x 

X = 2P = v(l-x), 
7 = vx. 



(4.23a) 

(4.23b) 
(4.23c) 



We find that our numerically determined exponents for 
< s < 1 obey these relations to high accuracy, as 
demonstrated in Fig. 1161 which compares numerical ob- 
tained exponents /3, 7 and 1/8 to values predicted by 
hyperscaling relations. Numerical values for representa- 
tive cases are given in Table 1 of Ref. 46. It follows from 
the relations Eqs. (4.23) and the asymptotic results Eqs. 
gSl and gTJ that for s -> 0+, 1/5 ~ 1 - 2s, A ~ 1/s 
and 7 ~ 1, while for s — » 1 _ , 1/S ~ i(l — s), A cx — s 
and 7 oc 1/y/l — s. 

We now turn briefly to the regime — 1 < s < through- 
out which the system is in the localized phase for any 
g > 0. Figure El shows the static local susceptibility 
Xioc(T;uj — 0) versus T for bath exponents s = and 
s = —0.8. The observed behavior is typical of any s in 
the range — 1 < < 0. As indicated by the flow dia- 
gram Fig. Hfb), the g — (Kondo) spectrum is followed 
down to progressively lower temperatures as g — > before 
crossing over directly to localized Curie-like behavior, Eq. 
iglTTJl. Consequently, for <C T K , xi oc (T;uj = 0) obeys 
the simple scaling form 



T K Xloc(T;U = 0) = <f) StL 7fT 



(4.24) 



(b) s = -0.8 
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FIG. 17: Static local susceptibility xioc(T;u = 0) versus T 
for J = 0.5 and (a) s — and (b) s = —0.8. The legends refer 
to the bosonic coupling g. Qualitatively similar behavior is 
observed in either case: as g — > there is a direct crossover 
from g — Kondo behavior for T > T,, to localized Curie- 
like behavior for T<T„ The localized-phase crossover scale 
T* vanishes as g —* — see Fig. 1181 



with no TVdependent prefactor, and v = X. Figure IT51 
shows Mi oc (0) 2 as a function of g 2 for the s values spec- 
ified. In all cases we find that Mi oc (0) 2 vanishes contin- 
uously with 



2\A 



M loc (0Y cx (g') 



(4.25) 



The inset plots the inverse of the exponent A versus bath 
exponent s. We find A = v ~ 1/s as s — * 0~, as found 
for s -> 0+ (see Fig. lTU|) . 



2. Dynamical susceptibility 

The local dynamical susceptibility xioc(w) can be com- 
puted using the NRG via its imaginary part 



with x ;; c » - xL(-«)> B « 



-El; 



(4.26) 
Em, and Z = 



Xloc V 



1 being the partition function. For T = 0, 
H = ^^|(n|5 2 |0)| 2 <5(|c| -Ko) sgn( W ) (4.27) 



with an implicit summation over ground states |0) in 
cases where the degeneracy Zq > 1 . Because only a small 
fraction of states are kept from one NRG iteration to the 
next, it is iteration N that contains the most appropri- 
ate information for dynamics at energy scales DAr N / 2 . 
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FIG. 18: M loc (T = 0) 2 = lim T ^o TxiocC? 1 ; u = 0) versus 
bosonic coupling g 2 for the s < bath exponents specified in 
the legend. For — 1 < s < the system is always localized for 
any g > and M\ oc {T — 0) 2 is found to vanish continuously 
as g — > 0. For j 2 C 1 the data can be fit to the form M\ oc (T = 
0) 2 oc (g 2 ) X - The inset shows the exponent A versus s, which 
diverges as A ~ 1/s for s — > 0~. See text for discussion. 



Thus Eq. (|4.26l) is evaluated at each NRG iteration for 
a frequency window spanning ujn — aDA~ N / 2 . The re- 
sult is a discrete set of delta-function peaks that must be 
broadened to recover a continuous spectrum. We follow 
well-established procedure^ applying Gaussian broad- 
ening on a logarithmic scale, viz 



5(uj - uj n ) 



exp 



(In a; — lna;„) z 



b 2 



(4.28) 



taking the broadening parameter 0.3 < b < 0.5. The 
static local susceptibility (considered in Sec. IIVD II) fol- 
lows via Hilbert transformation as 



Xioc(T;uj = 0) 



Xloc( £ > r ) 

7r e 



(4.29) 



Zero temperature. — Figure ED shows x(oc( w ) at T = 
on a logarithmic scale for a sequence of increasing cou- 
plings g < g c . For zero coupling to the bosonic bath, 
the familiar Kondo-model result obtains: Xk> c ( w ) ^ w 
for a; <C Tr-. However, for any nonzero g in the Kondo 
phase, x(o C (w) oc uj s as cj — * 0, corresponding to power- 
law relaxation xioc(i) oc i~( 1+s ) for long times. 

At the QCP (g — g c ) we find the scale-free form 



TkxL(^ 9 = 9c T = 0) cx C^j " Bgn(w), (4 



30) 



for w -C Tr - , with 




FIG. 19: Imaginary part of the local dynamical susceptibility 
X\ oc {u;T — 0) for bosonic bath exponent s = 0.5, Kondo 
coupling J = 0.5 and a range of g < g c ~ 1.216 in the Kondo- 
screened phase. The circles correspond to g = while the 
legend specifies g c — g for the remaining data. As g g^~ the 
curves follow the quantum-critical form for T* <C OJ <C Tr-. 
The dashed line is a fit to the latter [Eq. I|4.30|l ] which yields 
an exponent y — s to within less than 0.5%. See text for 
further comments. 



obeyed for all < s < 1 to within our estimated numer- 
ical error, which is typically less than 1%. This is illus- 
trated in Fig.EOI where x\o C {T; w = 0) and x" oc {lu] T = 0) 
are plotted on a logarithmic scale for bath exponents 
s = 0.2 and s — 0.8, and g = g c in either case. The 
result x — y is reflected in the equality of the slopes of 
the fits (shown as dashed lines). 

For small deviations from g — g c the critical power 
law behavior is cut off below a crossover scale T* that 
vanishes at the QCP according to Eq. I|4.5|l . Thus, the 
correlation length exponent v can also be extracted from 
x(o C (w;T = 0) and agrees with the results of Fig. ITUlto 
within our estimated numerical error. In the vicinity of 
the QCP we find the zero-temperature scaling form 



TkxI' oc (u;T = o) 



(4.32) 



y = x = s 



(4.31) 



for T* <C Tk, where p — K or L denotes the phase. For 
z ^> 1, 4> s ,p{z) = z~ s for either phase and thus Xi oc ( w ) 
is independent of T* for u»T, (see, e.g., Fig. ^J), fol- 
lowing the quantum-critical form TkXIoc(^) oh (lo/Tk)~ s 
before ultimately being cut off by the nonuniversal scale 
Tk/T*. Such scaling behavior is clearly illustrated in Fig. 
1211 which shows scaling for g — > g~ . 

Finite temperature. — The result y = x implies that 
the relaxation rate is linear in temperature, characteris- 
tic of critical physics below the upper critical dimension. 
Quantum-critical fluctuations are cut off on a timescale 
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FIG. 20: Critical X "oc(^9 = ffc,T = 0) versus w [xi oc (T;3 = 
ffcjW = 0) versus T] for J = 0.5 and s = 0.2 (filled symbols) 
and s = 0.5 (open symbols). It is seen clearly from the equal- 
ity of the slopes that the exponent y — x. Dashed lines are 
fits to the data and yield exponents y = x — s to within the 
estimated numerical uncertainty. 
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FIG. 22: Imaginary part of the dynamical local susceptibility 
xIo C (w,T) versus u), at the critical coupling g — g c for s = 0.5 
and a sequence of increasing temperatures T/Tk ~ (dotted 
line), 10 -5 , 10 -4 and 10 -3 . Temperature is the single low- 
energy scale, setting the location of the peak maximum. The 
inset shows data on a logarithmic scale for T/Tk — 10~ 2 , 
10" 4 , 10" 6 and 10" 8 . The T = data are show n as a dashed 
line and follow the critical power law Eq. (I4.30H with y = 0.5. 




FIG. 21: (Color online) Scaling of the T = dynamical susc- 
ceptibility for s=0.5 and g — > — see Eq. 14.3211 . The data 
were obtained for J = 0.5. The legend specifies g c — g and T» 
is here defined as the location of the peak maximum in Fig. 1191 
marking the crossover from quantum-critical uj~ s behavior for 
w > T, to io a for w < T,. For the data presented, T* varies 
over ten orders of magnitude. Excellent data collapse is ob- 
served up to scales set by the bare Kondo temperature Tk /T* . 
The dashed line is a simple fit taking 4> s> k{z) oc z s /(l + z 2s ), 
which does remarkably well at capturing the entire scaling 
curve. 



r = h/kgT. Thus, for nonzero T <C Tr- the quantum- 
critical form Eq. I|4.30|l is cut off for frequencies to « T, as 
is shown clearly in the inset to Fig.[221 For uj/T < T K /T 
this leads to ui/T scaling of the form 



TxxLc^^T^g = 9c) 



T 



0s (I) • (4.33) 



with 4> s (z) = z~ s for z » 1. The numerically obtained 
QCP scaling curve is shown in Fig. [53] as the thick line. 
Though the NRG method is unreliable for ui < T (Refs. 

andl57t) due to truncation error, we find a consistent 
description of the finite-T scaling properties of the model 
as discussed further below. See also Fig. 3 of Ref. l4fl 
which demonstrates u/T scaling for s — 0.2 and different 
values of T K . 

For any finite T», the local dynamical susceptibility 
obeys the full scaling form to be expected in the vicinity 
of an interacting QCP 19 : 



s.p 



Figure [23 plots ui/T scaling curves obtained via NRG 
calculations for several values of the ratio T*/T in the 
Kondo-screened phase of the model for s = 0.5. For 
lu/T > T« /r the curves follow the QCP (T» = 0) form 
[Eq. H4.33|) ]. which is of course independent of phase, 
p = K or L. Clearly, as T*/T — > the curves follow 
the QCP result more closely, i.e., 4> s (z) = ^ StP (z, 0). 
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E. Spectral function 
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FIG. 23: (Color online) uj/T scaling of the dynamical local 
susceptibility in the Kondo phase for s = 0.5 and values of 
T»/T specified in the legend. The QCP spectrum obtains 
for T,/T = 0. The tails of the curves fall onto the QCP 
form for T»/T < lj/T < T K /T, as shown clearly by the 
T„/T — 10 6 data. Here, T, is defined as the crossover scale 
in x{ qc (l>j;T = 0) from to~ s to lo s behavior. 



For the opposite extreme T*/T — * oo, the T = re- 
sult obtains when suitably rescaled (see Fig. I2T)) . such 
that & StP (z',z) — z~ s (f) s ^ p (z' j ' z) for z,z' — > oo. Since all 
data were obtained for finite Tk, the tail behavior is ul- 
timately cut off when lu/T increases to scales of order 
Tk/T, as seen in the figure. The T*/T dependence of 
the lj/T scaling spectrum has also been discussed in the 
context of the Kondo model with a pseudogap host den- 
sity of states p{e) cx |e| r , though the QCP of the latter 
model belongs to a different universality classics 

To summarize this subsection, we have studied in de- 
tail the local magnetic response for the Ising-anisotropic 
BFKM, extracting a range of critical exponents. For 
< s < 1 we find that the exponents satisfy hyperscal- 
ing relations and w/T-scaling characteristic of an inter- 
acting critical fixed point. All exponents associated with 
the local magnetic susceptibility can be determined by 
knowledge of just x and v. Furthermore, we find excel- 
lent agreement with corresponding results for the spin- 
boson model, confirming that the two models belon g to 
the same universality class. As pointed out in Ref. |4£| 
these results demonstrate the failure of a mapping 60,61 
of the SBM to a one-dimensional classical Ising model 
with long-range interactions falling off as l/r 1+s . The 
latter model is above its upper critical dimension for 
< s < 5 and the exponents there take mean-field val- 
ues with hyperscaling relations violated. However, in the 
range \ < s < 1 the SBM and classical long-range Ising 
model appear to be equivalent. See Refs.^g and fl9l for 
further discussion of this issue. 



We now turn to a discussion of the local single-particle 
spectrum, paying particular attention to the destruction 
of the Kondo resonance at the quantum-critical point. In 
this section we calculate the impurity spectral function 

A(w) = \ ]T \{n\di\m)\H{u - E n>m ) (e^™ + e~^") 

n.rn 

(4.35) 

for the Ising-symmetry Bose-Fermi Anderson model de- 
scribed by 

Him = y^(e d + \Uh d -<T)hda + Vy^X d 1 c k<r + H.c.) 

o k.a 

+ \g{fid ] -h di )Y,{^i + ^- q )- (4.36) 
i 

Here, ed is the impurity level energy, U > is the on- 
site Coulomb interaction, and V is the hybridization be- 
tween impurity and conduction band. In the limit where 
charge fluctuations at the impurity site are supressed, 
the Anderson model maps onto the Kondo model un- 
der a Schrieffer- Wolff transformation^ For particle- hole- 
symmetric impurity parameters = —\U, the mapping 
takes STq/ttU — + po^o with Tq = ttV 2 pq the hybridiza- 
tion width. 

The NRG approach described above, with the standard 
broadening procedure Eq. 14.28(1 applied to Eq. 14.35(1 . 
also leads also to a satisfactory description of A(uS) in 
the Bose-Fermi case, where the quantum phase transi- 
tion is manifest as a collapse of the central Kondo reso- 
nance. All results for A{lj) were obtained for NRG dis- 
cretization parameter A = 3, retaining N s = 800 states 
from one iteration to the next, and taking bosonic trun- 
cation parameter Nb = 8. Again, the half-bandwidth D 
is taken to be the energy unit unless stated otherwise. 
We use broadening parameter b = 0.3, which gives the 
most accurate recovery of the standard Bq = (K g ) 2 = 
Fermi-liquid-theory result A(ui = 0) = l/7rr (see, e.g., 
Ref. Ilfilh which is typically satisfied to within 2% within 
NRGS We do not vary b for Bq > 0. 

Fig. l2~4l compares the single-particle spectrum A(u>) - 
obtained for impurity parameters U — —26^ = 0.005 
[such that A(lj) = A(—u>)] and hybridization Tq = 
5 x 10~ 4 — for Bq = and two couplings < Bq < 
B Q , C ~ 0.047043(1). The coupling to the bosonic bath 
has two main effects. First, the Hubbard satellite bands 
are displaced to higher energies from their Bq = po- 
sitions, which in the Kondo limit (J< 1) are given by 
wh — ±5^- The shift is linear in Bq for small Bq, as 
is demonstrated in Fig. [53 Second, while the central 
Kondo resonance remains pinned to the Bq = Fermi- 
liquid-theory result A(uj = 0) = l/7rFo, which remains 
satisfied to within 2%, its halfwidth is significantly re- 
duced, as can be seen in the inset to Fig. |2] 

At the critical coupling Bq = B c , we find that the 
pinning of A(0) persists and the resonance has nonzero 
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FIG. 24: (Color online). Single-particle spectrum A(u) versus 
lu for s = 0.8, U = ~2e d = 0.005, r = 5 x 10~ 4 and Kondo 
phase bosonic couplings Bo < Bo, c — 0.04704 specified in the 
legend. The inset shows the central resonance on an expanded 
scale. See text for discussion. 



FIG. 26: Single-particle spectrum vs ui for bath exponent s = 
0.8, impurity parameters U = — 2ed = 0.005, hybridization 
strength To = 5 x 10 -4 , and discretization A = 3. Shown for 
Bo — Bo, c = (dashed line), 4 x ±10 -5 on (a) the scale set by 
the Kondo resonance; and (b) the scale of the on-site Coulomb 
repulsion U. The critical coupling Bo, c — 4.704 x 10~ 2 . 
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FIG. 25: Location of Hubbard band maximum as a func- 
tion of bosonic bath coupling Bo for s = 0.8, Fo = 5 x 10~ 4 , 
and the U = — 2ed values shown. The dashed lines are linear 
fits to the data for small Bo- 



width. This is illustrated in Fig. 1261 which shows three 
spectra: one at the QCP and one close to it in cither 
phase. Figure I26f a) plots the spectra on the scale of 
the central feature and shows clearly the pinning of the 
critical spectrum (solid line). For Bq > Bq, c the Kondo 
resonance collapses leaving a double peak structure with 
characteristic scale T* that vanishes as Bo — > Bq c . 
However, contrary to expectations, we find that A(0) 
is nonzero and decreases with increasing Bq. For fixed 
Bq > i?o, c we find that A(0) does not significantly vary 
with Nf, in the range 8 < Nf, < 16. Though the three 
spectra in Fig. 1261 show very different behaviors at low 
frequencies, on the scale of the Hubbard satellite bands 



[see Fig. I26f b)] they are indistinguishable. 

In the Kondo phase, the scale T„, which vanishes as 
Bq — > Bq c , is manifest in A(oj) as the crossover scale to 
the Fermi liquid form A(0) — A(uj) oc uj 2 as u> — > 0. 

Figure 1771 shows A(u>) for Bq — (dotted line) and a 
set of Bq values close to the critical coupling; A(u>) for 
Bq = Bq c is shown as a thick solid line. The crossover 
scale T* , at which the Fermi liquid "shoulder" is observed 
vanishes as Bq — > Bq c . On the localized side, the low 
frequency peak vanishes as Bq — > B^ c . Thus, regardless 
of phase, the critical spectrum is followed to progressively 
lower frequencies as the transition is approached. The 
correlation length exponent v(s) can be calculated from 
the crossovers in A(lo) described, and agrees well with 
the values determined for the BFKM in Sec. I1VI 

In the quantum-critical region \uj\ 3> T*, A(uj) is of 
form 



A(u = 0) - A{uj) cx \lo\ c 



(4.37) 



with A(oj = 0) = 1/7tFo satisfied for all < s < 1 to 
within 2%. To extract the power a we fit the critical 
spectrum (Bq = Bq^ c ) to the form 



A{uj) = [nT (l + k\u\ a )Y 



(4.38) 



shown as the thick dashed line in the inset to Fig. \77\ 
for s — 0.8. The fit is followed well over many orders 
of magnitude and is ultimately cut off by the Hubbard 
satellite bands. Though there remains a small numerical 
uncertainty in the fits, largely due to the fact that the 
zero- frequency pinning is only satisfied to within 2%, our 
results suggest that a = 1 — s for all < s < 1 with k a 
constant depending on s and Tk- 
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FIG. 27: (Color online) Single-particle spectrum A(uj) versus 
ui, shown on a logarithmic scale, for s — 0.8, U — —2t c i — 
0.005, To = 5 x 10~ 4 , and bosonic couplings Bo = (dotted 
line), B ,c (thick solid line), B , c ± 10" 5 and B , c ± 10~ 3 . 
See text for discussion. The inset shows a fit to the QCP 
spectrum of form A(lo) = [tiTo(1 + fcM")] -1 with a = 0.2, 
which is followed well over several orders of magnitude until 
it is cut off by the Hubbard satellite bands. 



Surprisingly, the key features described above are all 
present in the noninteracting limit (U = = 0) of the 
Anderson model, Eq. (|4.36|) . Figure |2"S1 shows the spec- 
tral function A(lo) obtained for s = 0.8 with (7 = 0, fixed 
ro = 5 x 10~ 4 , and a sequence of increasing bosonic cou- 
plings B . For B = (dotted line) the model is exactly 
solvable and it is well known that A(u>) has a simple 
Lorentzian form A(lo) — (Tq/tt)/ (lo 2 + Tq), which is re- 
covered well by the NRG broadening procedure. With 
increasing Bo < Bo :C w 0.0648 the central peak nar- 
rows rapidly, developing symmetric shoulders that sepa- 
rate into distinct bands redolent of the Hubbard satellites 
for the regular symmetric Anderson model (Bq = 0) for 
U = —2ed ^> 1. As found for U > 0, the position of 
the satellite bands increases linearly in Bo for small Bq, 
as shown in Fig. [23 At the critical coupling the cen- 
tral resonance has nonzero width and remains pinned at 
A(0) = l/7rFo. It is again possible to fit the critical spec- 
trum to the form Eq. (|4.38l) and, as for U = 0.005 (Fig. 
12711 . we find a — 0.2 for s — 0.8. This is demonstrated in 
the inset to Fig. [2H1 which shows data for Bo = 0.0648 
(Kondo phase) and 0.0650 (localized phase). In the fre- 
quency range shown the two data sets coincide — a reflec- 
tion of the fact that \u>\ 3> — and are well described 
by the fit discussed above, which is plotted as a dashed 
line. 

At the quantum-critical point itself we find that spec- 
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FIG. 28: (Color online) Single-particle spectrum A(lo) versus 
u> for non-interacting limit [7 = with To = 5 x 10 -4 , s = 0.8, 
and the Bo values shown in the legend. With increasing Bo < 
Bo, c ~ 0.0648 the central Lorentzian peak narrows rapidly, 
developing shoulders that separate from the parent peak to 
form distinct bands. In the quantum-critical regime (\uj\ S> 
T*) we find the same frequency dependence as for U > 0; the 
dashed line in the inset shows a fit of form Eq. 14.381 with 
a = 0.2 to data for B = 0.0648 (K phase) and 0.0650 (L 
phase). As for U > the localized phase spectrum contains 
a double peak structure about uj = with A(cu = 0) > 0. See 
text for further comment. 

tra scale according to 

ttT A(u;) = ^(uj/Bf 7 ), (4.39) 

as demonstrated in Fig. In the vicinity of the 

quantum-critical point single-particle spectra scale ac- 
cording to 

nT A(cj) = A, P Qr, ^) (4-40) 

where p = K or L denotes the phase. For the QCP, where 
T* = 0, and £?o, c oc TiT 6 , the form reduces to that of Eq. 
(|4.39|) . which scaling is demonstrated explicitly in Fig. 



V. RESULTS: OHMIC DISSIPATION, s = 1 

For any fixed Tk, the line of continuous quantum phase 
transitions in the g-s plane terminates at a Kosterlitz- 
Thouless-like transition at s = 1. As is to be ex- 
pected from the well-known mapping discussed in Sec. 
IIIBI the key physics of the BFKM with Ohmic dissipa- 
tion is already present in the anisotropic Kondo model. 
The Ohmic BFKM corresponds to an Ohmic spin-boson 
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FIG. 30: Schematic of NRG flow of effective couplings for the 
BFKM with Ohmic dissipation. The dashed line represents 
the phase boundary in the J-g plane. A crossover scale T* 
vanishes upon approaching the transition from the Kondo- 
screened phase. See text for discussion. 



FIG. 29: (Color online) Critical single-particle spectrum A{ui) 
for s = 0.8, U = 0.005 and 0.0025, e d = -\U, V = 5 x 10~ 4 , 
and Bo = _Bo, c in either case. The main figure shows scaling 
in terms of Sq^ 1- ^ oc Tr, while the inset shows the two data 
sets on an absolute scale. 

model with dissipation strength a containing contribu- 
tions from the Kondo term in the Hamiltonian, Eq. (|2.17|) 
with J± = J z = J , and from B , Eq. (|2.14() . For B = 
and < PoJq *C 1, the resulting SBM lies in the delo- 
calized phase with a = 1~; an exponentially-small renor- 
malized tunelling amplitude, corresponding to the Kondo 
scale, sets the scale for the low-energy physics. By in- 
creasing Bq > 0, the sytem may be further driven to- 
wards a Kosterlitz-Thouless transition at a = 1, where 
the renormalized tunelling amplitude vanishes, and even- 
tually into the localized phase a > 1. This is, of course, 
the essential physical behavior of the anisotropic Kondo 
model. In other words, the principal effect of coupling 
to an Ohmic bosonic bath is to introduce anisotropy into 
the corresponding Kondo model: PqJ± — PqJ z oc B . 

In this section we present results of direct calculations 
for the Bose-Fermi Kondo and Anderson models obtained 
using the NRG approach described in Sec. IIIII We begin 
by discussing the Kosterlitz-Thouless nature of the NRG 
flow of effective couplings. We then cast results for the 
static local susceptibility in terms of a noisy quantum 
box^ showing that perfect Coulomb-staircase behavior 
obtains above a critical coupling to the Ohmic bath. We 
conclude by discussing the single-particle spectral func- 
tion for the Bose-Fermi Anderson model with Ohmic dis- 
sipation, along with its universal scaling properties. 

For s — > 1~ the critical fixed point merges with the 
weak-coupling fixed point located at J = g = 0, as found 
in Ref. UX, A schematic NRG flow diagram for s = 1 is 
given in Fig. I3UI The phase boundary between Kondo- 
screened and localized phases is shown as a dashed line. 
Figure mi shows the numerically obtained phase bound- 
ary for the Ohmic case. As expected by analogy to the 
AKM, g\ oc J for J < 1. For g -> g~ the NRG level flow 
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FIG. 31: Phase boundary between Kondo (g < g c ) and lo- 
calized phases (g > g c ) as a function of Kondo exchange 
interaction J for Ohmic dissipation. For J < 1, j c ! a J 
(dashed line) . The inset shows the vanishing of the crossover 
scale T», with the latter obtained from NRG levels flows, as 
the phase boundary is approached from the Kondo phase for 
J — 0.05. The dashed line shows a fit to the form Eq. 15.1|l 
with ~ 0.958. 



is first towards the localized fixed point at J = g = 0. At 
a low-energy scale T* (large NRG iteration number N+), 
the NRG levels cross over to their stable Kondo-fixed- 
point values for s = 1. In contrast to the sub-Ohmic 
case < s < 1 , there is no crossover scale observed when 
approaching the transition from the localized side, cor- 
responding to a line of stable fixed points being reached 
for g > g c . 

The above comments are well illustrated by Fig. 1321 
which shows xioc(T;lu = 0) for s = 1, J = 0.5 and a 
range of dissipation strengths g spanning g c ~ 2.92. The 
observed behavior should be compared to that of Fig. 
E3for s = 0.5, which is typical of the sub-Ohmic case 
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FIG. 32: (Color online) x\oc(T; w = 0) versus T for the Ohmic 
case s = 1, J = 0.5 and bosonic couplings shown; g c ~ 2.92. 
On the Kondo side, a clear crossover scale exists, equivalent 
to the renormalized Kondo scale 0)> anc ^ vanishes con- 

tinuously as g — » (?~ . No such scale is present for the localized 
phase. Data obtained for A = 9, Nb = 12, and N 3 = 800. See 
text for further comment. 
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FIG. 33: Behavior of the renormalized Kondo scale X^o C (T = 
0, uj = 0) (oc T» for s = 1) and the order parameter M\ oc {T = 
0) for the Ohmic case s = 1 and J = 0.5. The dashed line is 
a fit of form Eq. 15. H . In the Ohmic case the order paramter 
undergoes a clear jump at the transition point g — g c ~ 2.92. 
See text for further discussion. 



s < 1. For Ohmic dissipation and g < g c , the impurity 
spin is quenched below a crossover scale T* that vanishes 
continuously as g — > g~ . The renormalized Kondo scale 
Xioc(0,0)~ 1 (oc T* for s = 1) likewise vanishes and we 



5 Z 



(5.1) 



Such a form follows from a poor-man's scaling analysis 
for the anisotropic Kondo modeljS which can be shown 
to yield 



Tk = D exp 1 



1 



\PqJz 



(5.2) 



which holds for <C 6 -C —PqJ z close to the transition 
with 5 2 = (poJ±) 2 — (poJz) 2 , and recognizing that g 2 — 
g 2 oc (poJ± — Po\Jz\) — S 2 /2J. A fit to the behavior Eq. 
(|5.1() is shown in Fig. [211 yielding a value for g c ~ 2.92 in 
good agreement with that obtained via NRG level flows. 
The corresponding fit for J = 0.05 is shown in the inset 
to Fig. 1^] The constants ao and ai vary as 1/ \fj for 
J <C 1, as expected from Eq. (|5.2|l . 

For T > T*, we find Xioc(T;w = 0) oc T" 1 with log- 
arithmic corrections. For g > g ci x\oc{T]uj — 0) oc 1/T 
shows Curie-like behavior characteristic of a freely fluc- 
tuating spin. However, as g — > g+, the order parameter 
Afioc(O) does not vanish continously, but rather jumps to 
zero at g = g c with the size of the jump being roughly 
90% of the full saturation value. 

The noisy quantum box. — The case of Ohmic dissipa- 
tion in the BFKM has recently been studied in the con- 
text of gate-voltage fluctuations in a quantum box. A 
quantum box consists of a small metallic grain (or dot) 
with a dense spectrum coupled to reservoir leads. Gate 
electrodes control the electrostatic energy of the dot, 



E c 



{Q-eNf 
2C 



(5.3) 



where N is proportial to the gate voltage, Q is the charge 
on the dot and C the capacitance. Degeneracy points at 
half-integer values of N (where Eq — Eq- c ) lead to a 
unit jump in the number of electrons on the grain at 
that point, and thus to the familiar Coulomb- blockade 
"staircase" behavior for (Q) versus N. However, the 
sharp charging steps are rounded due to quantum fluctu- 
ations, as shown theoretically by Matvee\i§i by mapping 
the problem to a two-channel Kondo model. Le Hur has 
recently shown^ that the supression of quantum charge 
fluctuations due to the dissipative effects of a fluctuating 
gate voltage can restore perfect Coulomb-staircase be- 
havior. This has been discussed in the context of a Bose- 
Fermi Kondo model, where the "spin" states of the im- 
purity correspond to the two charge states (orbital pseu- 
dospin) of the grain close to a degeneracy point. Con- 
sidering charging states Q — and Q = 1 leads to the 
identification (S z ) = (Q) — |. The electrons themselves 
are taken to be spinless, thus assuming the presence of a 
strong in-plane magnetic field. The gate voltage, which 
controls deviations away from the degeneracy point, i.e., 
differences in the energies of the spin states, plays the 
role of the local magnetic field h in the impurity problem. 
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FIG. 34: (Color online) Average charge on dot {Q) versus 
gate voltage N, for the noisy quantum dot system with (a) 
Ohmic and (b) sub-Ohmic dissipation. The legends specify 
the dissipation strength g in the BFKM. For weak dissipation 
the curves (Q) versus N are linear through the degeneracy 
point N = | due to Kondo screening of local moments. In 
the Ohmic case (a), the curves undergo a large jump at N = | 
for g > g c ~ 2.92. For g — g c ~ 0.566 in the sub-Ohmic case 
(b), (Q) — | vanishes in a power-law fashion at N = i, while 
for increasing g > g c it exhibits a jump of a magnitude that 
increases continuously from zero. 



The Kondo term represents tunneling from lead to grain 
and vice versa, while the bosonic bath, which is Ohmic, 
represents the voltage fluctuations. It then follows that 



d(Q) d(S z ) 
oc 



dN 



i:)h 



(5.4) 



The static susceptibility xioc(0,0) thus measures the 
slope of (Q) versus N at the degeneracy point N — 
■i. Our results indeed show a perfect restoration of 
the Coulomb staircase for large enough coupling to the 
bosonic environment. Figure 1547 a) shows (Q) versus N 
plots for the Ohmic case where perfect Coulomb staircase 
behavior is restored for g > g c . That there is a jump in 
(Q) for g — g c itself stems from the Kosterlitz-Thouless 
nature of the transition discussed above. For sub-Ohmic 
baths, by contrast, a continuous critical curve obtains at 
g = g c , with a slope that is infinite at N = ^. For increas- 
ing g > g c , the size of the jump in (Q) at N — \ increases 
from zero, consistent with the behavior of M(0) in Fig. 
1551 Our NRG results, which solve the Ohmic BFKM di- 
rectly, are consistent with the picture described by per- 
turbative RG£ and by a bosonic NRG treatment following 
the mapping of the BFKM to a spin-boson models 

We now turn to a discussion of spectral functions for 
the Bose-Fermi Anderson model with Ohmic dissipation. 
Figure 1551 shows spectra for fixed U = —2ed = 0.005 and 
To = 5 x 10~ 4 , for which the corresponding J ~ 0.25, 
and a sequence of increasing dissipation strengths Bq. 
There are several key differences compared to the sub- 
Ohmic case. First, and most obviously, the spectrum at 
the critical coupling Bq = Bq^ c ~ 0.33 is not pinned to 
the Kondo phase value A(u> = 0) = (7rFo) _1 , but is much 



FIG. 35: (Color online) Single particle spectrum A(uj) versus 
lj for the Ohmic case s — 1 and U = — 2ed = 0.005, To = 
5 x 10 -4 , and the So values shown. See text for comments. 



reduced. Furthermore, there is essentially no frequency 
dependence to the critical spectral function below scales 
on the order of the Hubbard satellite bands. For the 
parameters considered in Fig. [5J3 A(uo — 0) ~ 0.5 at 
the critical coupling. This form is followed by the Bq = 
0.32 data down to T* ~ 10" 120 which sets the scale for 
crossover to Kondo-phase behavior. 

Second, the vanishing of the crossover scale T», de- 
scribing, e.g., the vanishing of the Kondo resonance half- 
width, is well fit by 



InT* oc ao 



a i 



V Bq, c - Bq 



(5.5) 



as expected from the result Eq. I|5.1|l for the Ohmic Bose- 
Fermi Kondo model. 

Third, the Kondo-phase spectra obey a simple scaling 
form 



ttT A(uj) = <J) S 



-l.K 



(5.6) 



as shown in Fig. I36f b) . The spectra there plotted, which 
have vastly different crossover scales (see, e.g., Fig. 13511 . 
collapse to the single universal curve Eq. I|5.6|l . which 
form is followed up to the nonuniversal energy scale set 
by the Hubbard satellite bands, ojh/T*- For the data 
shown, excellent data-collapse is seen up to uj/T* ~ 10 7 . 

We note that calculations of the T-matrix in Ref.0, dis- 
cussed in the context of the quantum box, are in qualti- 
tative agreement with the above. 



VI. SUMMARY 

We have described an extension of Wilson's numeri- 
cal renormalization group method appropriate for Bose- 
Fermi quantum impurity models. We have focused in 
particular on the Ising anisotropic Bose-Fermi Kondo 
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36: (Color online) (a) Vanishing of crossover scale T* for 
-> Bg c , according to Eq. 15.51 (dashed line), (b) Scaling 
of the single-particle spectrum nToA(uj) for the dissipation 
strengths Bo shown in the inset. Though the values of T* 
in this data set vary across ten orders of magnitude, a single 
universal scaling curve [Eq. 15.61 is followed up to the nonuni- 
versal scale set by the Hubbard satellite bands. 



model. The approach gives an excellent description of 
the critical properties of this model, which is of relevance 
to certain dissipative mesoscopic qubit devices and to a 
treatment of the anisotropic Kondo lattice model within 
an extension of dynamical mean-field theory^ For sub- 
Ohmic bosonic bath exponents < s < 1 we find a 
continuous quantum phase transition between Kondo- 



screened and localized phases. The transition is governed 
by an interacting quantum-critical point, leading to hy- 
perscaling of critical exponents and w/T-scaling in the 
dynamics. We have shown directly that the Kondo reso- 
nance in the spectral function is destroyed at the critical 
point. The case of Ohmic dissipation, which leads to a 
Kosterlitz-Thouless like transition, was also discussed. 

While we have specifically studied the case of Ising 
symmetry bosonic couplings, the approach is not per se 
limited to that case. It is straightfoward, albeit com- 
putationally more demanding, to add additional bosonic 
bath(s) as appropriate for the xy-anisotropic or isotropic 
Bose-Fermi Kondo models. These cases will be discussed 
elsewhere. Furthermore, our approach could be inte- 
grated with the recently developed time-dependent NRG 
methociSiS& to enable the study of nonequilibrium dy- 
namics for a wider range of quantum impurity systems. 
Other models that might be tackled using the present ap- 
proach include an Anderson-Holstein-type model (where 
the impurity charge density couples to a bosonic bath) 
relevant to the study of electron-phonon interactions in 
strongly correlated electron systems, and novel Bose- 
Fermi impurity problems that may arise in the treatment 
of ultracold mixtures of bosonic and fermionic atoms 
trapped in optical lattices. These applications are left 
for future study. 
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High Performance Computing Cluster; we thank C. Tay- 
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